Experimental Signatures of Critically Balanced Turbulence in MAST 
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Beam Emission Spectroscopy (BES) measurements of ion-scale density fluctuations in the MAST 
tokamak are used to show that the turbulence correlation time, the drift time associated with ion 
temperature or density gradients, the particle (ion) streaming time along the magnetic field and 
the magnetic drift time are consistently comparable, suggesting a "critically balanced" turbulence 
determined by the local equilibrium. The resulting scalings of the poloidal and radial correlation 
lengths are derived and tested. The nonlinear time inferred from the density fluctuations is longer 
than the other times; its ratio to the correlation time scales as v~°' where v»i = ion collision 

rate/streaming rate. This is consistent with turbulent decorrelation being controlled by a zonal 
component, invisible to the BES, with an amplitude exceeding the drift waves' by ~ v~ t ' 8 . 



Introduction. Microscalc turbulence hindering energy 
confinement in magnetically confined hot plasmas is 
driven by gradients of equilibrium quantities such as 
temperature and density. These gradients give rise to 
instabilities that inject energy into plasma fluctuations 
("drift waves") at scales just above the ion Larmor scale. 
The most effective of these is believed to be the ion- 
temperature-gradient (ITG) instability A turbu- 
lent state ensues, giving rise to "anomalous transport" of 
energy Q. It is of interest, both for practical considera- 
tions of improving confinement and for the fundamental 
understanding of multiscale plasma dynamics, what the 
structure of this turbulence is and how its amplitude, 
scale (s) and resulting transport depend on the equilib- 
rium parameters: ion and electron temperatures, density, 
angular velocity, magnetic geometry, etc. 

Fluctuations in a magnetized toroidal plasma are sub- 
ject to a number of distinct physical effects, which can be 
thought about in terms of various time scales such as the 
drift times associated with the temperature and density 
gradients, the particle streaming time along the magnetic 
field as it takes them around the torus toroidally and 
poloidally, the magnetic (VB and curvature) drift times 
of particles moving across the field, the nonlinear time 
of the fluctuations being advected across the field by the 
fluctuating E x B velocity, the time between collisions, 
the shear time associated with plasma rotation. Some 
of these time scales and, consequently, the correspond- 
ing physics may be irrelevant, while others play a cru- 
cial role for the saturation of the linearly unstable fluc- 
tuations. There has been agrowing understan ding 
driven largely by theory [6H9|, observations [10h12| and 
simulations of magnetohydrodynamic I3-15| and kinetic 
0, EH plasma turbulence in space, that if a medium can 



support parallel (to the magnetic field) propagation of 
waves (and/or particles) and nonlinear interactions in the 
perpendicular direction, the turbulence in such a medium 
would normally be "critically balanced," meaning that 
the characteristic time scales of propagation and nonlin- 
ear interaction would be comparable to each other and 
(therefore) to the correlation time of the fluctuations. 
This means that the turbulence is not weak and not two- 
dimensional, unless specially constrained to be so Q. 

Beam Emission Spectroscopy (BES) measurements of 
density fluctuations in tokamak plasmas 17H2l| have 
made it possible to probe ion-scale turbulence in these de- 
vices directly. In this Letter, we use such measurements 
in the MAST tokamak, along with the local equilibrium 
parameters calculated by other diagnostics, to estimate 
and compare the characteristic time scales of the tur- 
bulent fluctuations in the energy-containing range. We 
obtain, for the first time, direct evidence that the corre- 
lation, drift and parallel streaming time scales are indeed 
comparable across a range of equilibrium parameters (cf. 
22l |23| ) and that the magnetic drift time is part of this 
"grand critical balance" as well. We also find indirect 
evidence that the decorrelation rate of turbulence is con- 
trolled by a zonal component whose relative importance 
to the drift-wave-like fluctuations scales with the ion col- 
lisionality. 

Before presenting this evidence and its implications 
(e.g., dependence of the correlation lengths on equilib- 
rium parameters), let us describe how it was obtained. 

Experimental data and its analysis. During the 2011 
campaign, density fluctuation data from the BES diag- 
nostic |21[ on MAST were collected in a variety of dis- 
charges (including L- and H-modes and internal trans- 
port barriers). Here we report the data from 39 neutral- 
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FIG. 1. An example of the correlation function in the 
poloidal-temporal plane, C (Ax = 0, AZ, At). This data was 
taken at r — 30 cm, toroidal rotation speed was 11$ = 10 km/s 
and magnetic pitch angle a = 20° . The direction of maximum 
correlation is the direction of the magnetic field (dashed line) . 

beam heated "double-null-diverted" discharges, with no 
pellet injection and no resonant magnetic perturbations. 
The BES system on MAST collects photons from a 
2D array of 8 radial x 4 vertical locations in the out- 
board midplane of the tokamak, with 2 cm separation 
between the adjacent channels in either direction. The 
detected photon intensity (mean + fluctuating, / + SI) 
is used to infer, at each location, the density fluctuation 
level Sn/n = (1//3) (SI/I) [T3], where (3 depends on the 
mean density n and is estimated based on the Hutchin- 
son model [2J] (dependence on the mean temperature is 
weak). As the BES array was moved radially for differ- 
ent discharges, our database contains cases with radial 
viewing positions 10 cm < r < 50 cm from the magnetic 
axis (the minor radius of the plasma is w 60 cm). 

Local equilibrium parameters are measured by stan- 
dard diagnostics: mean electron densities n e and tem- 
peratures T e by the Thomson scattering system [25| , im- 
purity ion (C 6+ ) mean temperatures (assumed to equal 
the bulk ion temperature T{) and toroidal flow velocity 
U$ by the Charge eXchange Recombination Spectroscopy 
(CXRS) system [26|, local magnetic pitch angle a by the 
Motional Stark Effect (MSE) system [27|, and further 
equilibrium magnetic field information is obtained from 
pressure- and MSE-constrained EFIT equilibria pjj ]. 

We filter the BES data to the frequency interval 
[20, 100] kHz and calculate the spatio-temporal correla- 
tion function 

C (Ax, AZ, At) = 

(51 (x, Z, t) SI (x + Ax,Z + AZ, t + At)) 
y/{SP (x, Z, t)) (SP (x + Ax,Z + AZ, t + Al)) ' 

where x, Z and t are the radial, vertical and time coor- 
dinates, respectively, and Ax, AZ and At are the cor- 
responding channel separations and the time lag; (•) is 
the time average over 5 ms periods. At Ax = AZ = 
0, the auto-covariances (SI (x, Z, t) 81 (x, Z,t + At)) con- 
tain not only the physical signal but also photon and elec- 
tronic noise. We remove this effect by applying LED light 



to the BES channels, obtaining 150 different DC levels of 
BES signal from to 1.5 V, calculating the noise auto- 
covariance CV (At) at each DC level with the same band 
frequency filter of [20, 100] kHz, then finding C N (At) 
whose DC level of the signal matches the DC level of the 
BES data from the MAST discharges, and subtracting it 
from the calculated auto-covariances. From the correla- 
tion function (JlJ (illustrated in Fig. [IJ , we calculate the 
local characteristics of the density fluctuations. 

The fluctuation level at each radial location is ob- 
tained from the (noise-subtracted) auto-covariancc func- 
tion Sn/n = (1/(3) y/(SP(x,Z,t))/I at all 32 locations 
and then averaged over the four poloidally separated 
channels at the same radial location. 

The correlation length l y in the direction parallel 
to the flux surface and perpendicular to the mag- 
netic field is obtained from the vertical (poloidal) 
correlation length £z via £ y = izcosa, assuming 
that the parallel correlation length is sufficiently long: 
£\\ 3> l y tan a. The correlation length lz is esti- 
mated using four poloidal channels at each radial loca- 
tion (the top channel is the reference channel) by fit- 
ting C (Ax = 0, AZ, At = 0) to the function f z (AZ) = 
Pz + 0-- Pz) cos [2irAZ/£ z ] exp [- \AZ\ /£ z ], where p z 
is a fitting constant that serves to account for global 
structures such as coherent MHD modes (for which 
C (Ax = 0, AZ = oc, At = 0) = p z / 0). In choos- 
ing fz ('AZ), we assumed wave- like fluctuations in the 
poloidal direction [l8| (drift- wave turbulence), with the 
wavelength and correlation length comparable to each 
other. It is not possible to distinguish meaningfully be- 
tween the two with only four poloidal channels. Assum- 
ing wave-like structure is essential as in most cases, we 
find that C (Ax = 0, AZ, At = 0) goes negative and/or is 
non-monotonic over the vertical extent of the BES array. 

The radial correlation length £ x is estimated using 
eight radial channels at each poloidal location (the fourth 
channel from the inward side is the reference channel). 
The correlation function C (Ax, AZ = 0,At = 0) is fitted 
to the function f x (Ax) = p x + (1 - p x ) exp [- |Ax| /£ x ], 
where p x plays the same role as pz did for fz . The values 
of £ x from four poloidal locations are averaged, assuming 
that the radial correlations do not change significantly 
within the poloidal extent of the BES array. Because we 
have to use the entire array to estimate l x , the number 
of data points for £ x is 8 times smaller than for £ y . 

To estimate the correlation time t c , we use the fact that 
the fluctuating density patterns are advected poloidally 
past the BES array with an apparent velocity vbes = 
U,), tan a due to the toroidal rotation velocity U$ [29( . 
We fit C (Ax = 0, AZ, At = At pcak (AZ)) taken at the 
time delay At poa k (AZ) when the correlation function is 
maximum at a given AZ [30l ] , to the function f T ( AZ) = 
exp [— |A< pca k (AZ)| /t c ]. The reliability of this method 
relics on the temporal decorrelation dominating over 
the parallel spatial decorrelation, viz., we require r c <C 
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£ii cos a/Uij). Anticipating the critical balance assumption 
t c ~ £\\/v t hi 0, where v t hi = \J ^Ti I mi is the ion thermal 
speed, and denoting the Mach number Ma = U^/vthi, we 
estimate that the fractional error in r c is ~ Ma/ cos a, 
which was never more than 20% in the MAST discharges 
we used. 

The four quantities dn/n, £ y , £ x and t c are calcu- 
lated at 8 radial locations (except £ x ), every 5 ms for 
all 39 discharges. All the fits described above are ob- 
tained via the mpf it procedure (3l| . We consider a data 
point unreliable and remove it from the database if (i) 
/ < 0.3 V (the signal-to-noise ratio is too low); (ii) the 
estimated correlation lengths are smaller than the dis- 
tance between the channels, £ x or l y < 2 cm; (iii) the 
assumption that plasma rotation is mostly toroidal is 
suspect, viz., \(vbes — tana) /vbes\ > 0.2 (see ref. 
[29l| ) . where vbes is calculated at each radial location 
using the cross-correlation time delay (CCTD) method 
[30l |: (iv) the estimated error in the calculation of vbes 
is > 20%; (v) pz or p x > 0.5. The last two exclusion cri- 
teria pick out the cases when MHD modes are too strong; 
they are known to degrade the reliability of the BES data 
[29l | . The remaining database contains 448 points. 

Correlation time vs. drift time. The turbulence can 
be driven by radial gradients in the mean ion and elec- 
tron temperatures Ti^ e and density n. Denoting L^} = 
|VlnTj ie | and L^ 1 = |Vlnn|, the associated time scales 
are the inverse drift frequencies: 
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where p, ie = Vthi,e/^i,e are the ion (i) and electron (e) 
Larmor radii, «thi,e = y/^Ti^/m^e the thermal speeds 
and rii^ e = eB/rrii^c the Larmor frequencies. To es- 
timate the drift times, we need information about the 
local equilibrium (T^ e , Lt; e L„, B) and the correlation 
length £ y , calculated from the poloidal BES correlations. 

In Fig. EJa) , we compare the drift times with the cor- 
relation time t c calculated from the spatio-temporal BES 
correlations. We find that r* = (0.7 ± 0.3)r c , where 
r* = min{r*i,r*„} and the spread is calculated as the 
root mean square deviation from the mean value. The 
scaling holds over an order of magnitude in either time 
scale. Thus, the turbulence appears to be driven by the 
larger of the ion temperature or density gradient [Hf . We 
find no clear correlation of r* e with r c , or with any of the 
other time scales discussed below. 

Critical balance. The standard argument behind the 
critical balance conjecture is causality @: two distant 
points on a field line cannot stay correlated if information 
cannot be exchanged between them over a turbulence 
correlation time. Assuming information travels at Vthi, 
one gets £\\ ~ VthiT c . This cannot be checked directly 
because there are no diagnostics capable of measuring 
£» on MAST 33j. Considering that the inboard side of 
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FIG. 2. (a) Drift time r* = (£ y /pi)L*/v t hi vs. correlation 
time t c ; (b) streaming time r st = A/i>tM = (B / B p )nr / v t hi 
vs. t c ; (c) magnetic drift time r M = (£x/pi)R/v t hi vs. r c ; (d) 
perpendicular velocity shear time r s h = [(B p /B)dU r f,/dr]~ 1 
vs. r c . In all cases, the color of points represents r\i — Ln/L^ ■ 



the torus is a region of "good" (stabilizing) curvature, not 
much turbulence is expected there .so we assume that, at 
the energy injection scale, £\\ ~ A Q, where the distance 
along the field line that takes a particle from the outer to 
the inner side of the torus is A = irrB / B p (r is the minor 
radius at the BES position on the outer side and B p the 
poloidal component of the magnetic field) 34(. Then 
critical balance means that r c should be comparable to 
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the ion streaming time (the first two equalities are its 
definition, the last an assumption). Indeed, we find 
r st = (0.8 ± 0.3)t c (see Fig. 02b)). 

The balance r st ~ r* implies that the poloidal corre- 
lation scale is £ y /pi ~ A/L», where L* = minjXT^ , L n } 
This is tested in Fig.EJa), showing that while the two 
quantities are certainly of the same order, we do not have 
enough of a range of equilibrium parameters to state con- 
clusively that this theoretically predicted scaling works. 

Magnetic drift time and radial correlation scale. The 
time scale of the magnetic (VB and curvature) drifts is 
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where we have assumed that the scale length of the back- 
ground magnetic field is R (major radius at the viewing 
location) and £ x < £ y (this will shortly prove correct). It 
is clear that this scale cannot be shorter than r c because 
damping due to the drift resonance would eliminate such 
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FIG. 3. (a) Poloidal correlation length £ y /pi vs. A/L„; (b) 
Radial correlation length i x /pi vs. A./R. Color as in Fig. [2] 



FIG. 4. (a) The nonlinear time associated with density fluc- 
tuations, r^} z , vs. the correlation time r c ; (b) their ratio vs. 
normalized ion collision rate v*i = fuT^t- Color as in Fig. [2] 



fluctuations. While magnetic drift physics may matter 
(in a torus, curvature contributes to the ITG drive Q), 
it does not have to affect scalings, as, for example, it 
would not in a slab and as it did not in the numerical 
simulations of [H[ . In contrast, Fig. [^c) shows that in 
the MAST discharges we have analyzed, tm is not negli- 
gible and scales with r c , similarly to t* and r st [Hj]. We 
find tm = (1.6±0.7)t c . Thus, a "grand critical balance" 
appears to hold in MAST, viz., r c ~ t* ~ r st ~ tm- 

This suggests that the balance of all relevant 
timescales determines correlation scales of the turbulence 
in all three spatial directions. Indeed, balancing tm ~ T st , 
we find the radial correlation scale £ x /pi ~ A/JJ, the scal- 
ing tested in Fig. [3jb), with a degree of success. This 
means that the density fluctuations we are measuring in 
MAST are not isotropic in the perpendicular plane, but 
rather elongated in the poloidal direction £ y /£ x ~ RjL % 
(~ 5 in our data). Interestingly, this clashes with the 
reported approximate isotropy {l x ~ &y) both in Cyclone 
Base Case simulations [f| and in measured DIII-D turbu- 



lence (where 



1.4 [36[ and £ x does not appear to 



depend on B p [37|). Whether this is a difference between 
spherical and conventional tokamaks is not as yet clear. 

Nonlinear time. Since we know the fluctuation ampli- 
tude, we can directly estimate the time scale associated 
with the advection of the fluctuations (Su± ■ VSn) by the 
fluctuating E x B velocity Su±^ = cB x Vtp/B 2 . The 
electrostatic potential tp is not directly measured, but 
can be estimated assuming Boltzmann response of the 
electrons: Sn/n « eip/T e . This estimate ignores trapped 
particles and, more importantly as we are about to argue, 
also does not apply to ion-scale zonal flows (poloidally 
and toroidally symmetric perturbations of tp with 5n = 
*38l 3t|). Thus, the non-zonal nonlinear time is 
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Fig. Sta) shows that r^ z is always larger than r c (or the 
other time scales discussed above) and, furthermore, ob- 
served to have an inverse rather than direct correlation 
with it. Since turbulence clearly cannot be saturated by 
linear physics alone, this means that our estimate does 



not capture the correct nonlinear time. We conjecture 
that it is in fact the coupling to the zonal flows, invisible 
to BES (because their Sn = 0), that dominates over the 
nonlinear interaction between the drift-wave-like fluctu- 

It has long been 



ations represented by [38|, |4C 
suspected that the relative amplitude of the zonal flows 
compared to that of the drift waves depends on the ion 
collisionality [H We can test this expectation 

by assuming that r c is the characteristic time associated 
with the coupling of the drift waves to the zonal flows and 
so depends on their amplitude. Fig. |4jb) indeed shows 
a strong collisionality dependence: t c /t^ z ~ ^ 8±01 , 
where v„ = z/«T s t [50( (the ion collision time itself, u# , 
is at least an order of magnitude longer than the time 
scales that participate in the "grand critical balance"). 
If r^ 1 ~ (vthiPi/£ x £ y )e<p ZF /Ti, where tp ZF is the ampli- 
tude of the zonal potential, this result implies that the 



ratio of zonal to non-zonal component of the turbulence 
is ^ ZF /^ NZ ~ l /-°- 8±0 - 1 ' 



51] 



We note that this situation is qualitatively distinct 
from what is seen in numerical simulations of ITG tur- 
bulence far from the threshold where the drift-wave 
nonlinearity appears to dominate (t^ z ~ t c ). However, 
the turbulence in a real tokamak is likely to be close to 
marginal and so possibly in the state of reduced transport 
controlled by weakly-collisionally damped zonal flows [39| 
and usually associated with the so-called "Dimits up- 
shift" of the stiff-transport threshold 0, Q [H, H3] . 

Discussion. Our results support the notion that the 
statistics of turbulence are determined by the local equi- 
librium properties of the plasma. We find little correla- 
tion between the quantities reported above and the radial 
location [HiJ (note that we have limited our consideration 
to temporal and spatial scales and did not touch on the 
fluctuation amplitudes or transport properties, which do 
of course depend on radius). Our results also appeared 
insensitive to (i.e., not measurably correlated with) three 
other parameters that might in principle have proven im- 
portant: Ti/T e (varied between 0.5 and 2), the magnetic 
shear s = d\\\q/ dh\r (varied between —1 and 5) and the 
perpendicular component of the toroidal velocity shear 
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= (B p /B)dU t p/dr. In much of our data, r s h > r c ,T st 
(see Fig. Efd)), so it stands to reason that the statistics 
of the turbulence would not be dramatically affected; in 
the instances of r s h ~ T s t, the effect of r s h could not 
be isolated (55[. It would be interesting to investigate 
higher-rotation plasmas, as tT , when sufficientl y la rge, 



is expected to have a dramatic effect on transport [411.156- 
[6^ |: even in our database, there is in fact some evidence 
that velocity shear might raise the critical temperature 
gradients [631 ] . but we see no signature of this effect in 
the correlation properties of the turbulence. 

Conclusion. We have presented experimental results 
statistically consistent with a turbulent state in MAST 
set by the local equilibrium and in which the time scales 
of the linear drive, turbulence decorrelation, ion stream- 
ing and magnetic drifts are all similar and scale together 
as equilibrium parameters are varied. This "grand crit- 
ical balance" implies a three-dimensionally anisotropic 
turbulence, with parallel, poloidal and radial correla- 
tion lengths having different parameter dependences and 
£\\ 3> i y > £ x . Our results also suggest the presence of 
a zonal component with an amplitude ^~°' 8±01 greater 
than the drift-wave density fluctuations. 
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